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Interpreting Magnetic Variance Anisotropy Measurements in the Solar Wind 

J. M. TenBarge, 1 ^ J- J- Podesta, 2 and K. G. Klein and G. G. Howes 1 

^Department of Physics and Astronomy, University of Iowa, Iowa City, I A, 
USA 

2 ' Space Science Institute, Boulder, CO, USA 

The magnetic variance anisotropy (A m ) of the solar wind has been used widely as a method to identify the 
nature of solar wind turbulent fluctuations; however, a thorough discussion of the meaning and interpretation 
of the A m has not appeared in the literature. This paper explores the implications and limitations of using 
the Am as a method for constraining the solar wind fluctuation mode composition and presents a more 
informative method for interpreting spacecraft data. The paper also compares predictions of the A rn from 
linear theory to nonlinear turbulence simulations and solar wind measurements. In both cases, linear theory 
compares well and suggests the solar wind for the interval studied is dominantly Alfvenic in the inertial and 
dissipation ranges to scales kpi ~ 5. 
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I. INTRODUCTION 

In situ measurements of solar wind turbulence con- 
sistently show one dimensional magnetic energy spec- 
tra that obey a broken power law, typically having 
spectral indices —5/3 in the i nertial range and steep- 
ening in the dissip ation range (lAlexandrova et al.l [20 lit 
iBruno fc Carbond 120051: iTu fe Marschlll995l ). The iner- 
tial and dissipation ranges correspond to scales kr{ < 1 
and fcr,; > 1 respectively, where k is the wavenumber and 
Ti is the relevant ion kinetic scale, typically either the ion 
gyroradius or inertial length. 

Although measurements of the magnetic energy spec- 
trum are common, the spectrum alone does not provide 
direct insight into the nature of solar wind fluctuations. 
Since the solar wind turbulence is electromagnetic, the 
expectation is that the solar wind fluctuations will ex- 
hibit characteristics of the three basic electromagnetic 
plasma wave modes at the large scales of the inertial 
range: Alfven, fast m agnetosonic, and slow magnetosonic 
( Klein et al.l l2012bh . Similarly, the dissipation range is 
expected to be populated by the kinetic scale counter- 
parts of the three wave modes. 

Based upon a variety of metrics, the Alfven mode 
appears to b e the dominant wave m o de in the in- 
ertia! range (IBruno fe Carbonel |2005|: iHorburv et al 



20081: iHowes et al.l l2011at iPodestal l2009UTu fc Marsch 



19951 ). The composition of solar wind fluctuations 



in the dissipation range is less well constrained due 
to the dearth of high frequency measurements in the 
free solar wind necessary to probe this region, but 
recent observations suggest the kinetic Alfven wave 
( KAW) is the dominant mode in the free solar wind 
Bale et alj|2005l IChandran et alj|2005 iPodesta fc Garvl 



l201lUSahraoui et alJl2010USalem et al.ll2012h — the KAW 

is the dissipation range extension of the inertial range 
Alfven wave in the k± > ku region of wavenumber space, 



where parallel and perpendicular arc with respect to the 
local mean magnetic field, Bq . 



One of the commonly used metrics is the mag- 
netic variance anisotr opy (A m ), first introduced by 
iBelcher fc Davisl (|l97lh . The A m is a measure of the 
fluctuation anisotropy and has come to be defined as 



A m = (\o~B±\ 2 )/ (\5B\\ | 2 ), where angle brackets indicate 
averages, the SB are fluctuating quantities about the lo- 
cal mean magnetic field, and 8B\ is the total energy 
in the plane perpendicular to the local mean magnetic 
field. It is important to not confuse this quantity with 
the wavevector anisotropy inherent to and often discussed 
in plasma turbulence: the wavevector a nisotropy and the 
Am a re not directly related quantities ( Matthaeus et al.l 
1995). Physically, the A m can be viewed as a mea- 
sure of the magnetic compressibility of the plasma, Cii = 
(|<5E||| 2 )/(|<5B| 2 ), since A m + 1 = l/Cy. 



The version of the 

IBelcher fc Davisl (ll97ll) 



Ajr 

has 

i n recent pape r s usi n g the A 



first introduced by 
been expanded upon 



(iGarv fc Smith! 120091; lHamilton et al 



2012 



defin ed ab ove 
20081: 



Leamon et al.l Il998t IPodesta 



.He et al 
TenBared l2"oT ~ 



et al.l 
20121: 



a ' Electronic mail: jason-tcnbargc@uiowa.edu 



Salem et alj|2012l ISmith et alj|2006ll2012h . Weattempt 
here to provide a discussion that includes a theoretical 
basis for the interpretation of A m measurements in 
the solar wind. In JEH wc describe in detail the ex- 
pected behaviour of the A m of the three constituent 
electromagnetic wave modes in the solar wind inertial 
range and discuss their transition into the dissipation 
range. In £11111 we discuss the effect of superposing the 
three wave modes. ^IVI explores alternative procedures 
for constructing the A m , and compares the linear 
theory prediction from fjll] to fully nonlinear turbulence 
simulations as a means establish the validity of linear 
theory to nonlinear turbulence. <JV] reviews some of the 
recent uses of the A m to quantify the composition of 
solar wind fluctuations and presents new measurements 
of the A m from the Stereo A spacecraft. 
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FIG. 1: (Color online) Schematic diagram of wave mode 
transitions in wavenumber space. Black, red, and blue 
indicate Alfven, fast, and slow mode branches 
respectively. 



II. WAVE MODES 

We begin by enumerating the properties of the three 
linear wave modes which are the collisionless counter- 
parts to the Alfven and fas t and slow magne tosonic 
modes in compressible MHD ([Klein et al.ll2012b|) . Since 
the solar wind is a weakly collisional plasma, we focus 
here on the roots provided by the collisionless Vlasov- 
Maxwcll (VM) system of equations, which are a function 
of k\\,k±, ft, Ti/T e , and v t i/c, where v% = ITi/rrii is 
the ion (protons only) thermal speed. v t i/c = 1CP 4 and 
Tj = T e unless otherwise stated. Figure [T] presents a 
schematic diagram summarizing the nomenclature of the 
three wave branches in different regions of wavenumber 
space. 



A. Alfven Mode 



At scales kj_pi <c 1 and kudi 



< 1 in the free solar 



mode ( 


Bruno & Carbone 2005; 


Podesta 


20091 iPodesta & Garvl 


20101: iTu & Marsch 11995). where 



= Uti/Oj is the ion 
(proton) gyroradius, Oj = eBo/rriiC is the proton gy- 
rofrequency, di 



„2 , 



is the ion inertial 



length, and = iwriie /rrii is the ion plasma fre- 
quency. The turbulent energy cascade at these scales 
has been thoroughly explored in the literature, where 



the one dimensional perpendicular magnetic energy is 
observed to scale as Eb ± oc kj 01 with a wavevector 
anisotropy feii cx fc^. When the turbulence is in crit- 
ical balance — i.e, the nonlinear cascade rate is of or- 
der the linear frequency — the theoretical values for a 
and £ are expected to be a = 5/3 and £ = 2/3 
for th e model of Goldreich-Sridhar ( Goldreich fc Sridharl 
Il995l ) or a = 3/2 and £ = 1/2 for the dynamic align- 



ment model (|Boldvrevl 120051 120061) . Which model is 
correct is not completely settled. Solar wind obser- 
vations typically show a magnetic field spectrum with 
a ~ 5/3 and total (ki netic plus magnetic) e nergy spec- 
trum wit h a ~ 3/2 (iBoldvrev et all l201li IChen et al l 
2011alibl: iLuo fc Wul 120101 : IPodesta fc Borovskvl 120101: 



Wicks et all I2010T) . while MHD turbulence simulations 



sugg est the dynamic alignmen t model may be more cor- 
rect ijMason et al.ll2006l I2008T) . For the purpose of mod- 
elling the energy cascade, we assume the Goldreich- 
Sridhar model for simplicity. The choice of model does 
not significantly affect the behaviour of the A m - 

The properties of the MHD Alfven root are well un- 
derstood. However, the MHD solution of the Alfven root 
is incompressible since <5-B|| = 0, suggesting the Alfvcnic 
A m is unbounded. The full collisionless VM solution of 
the Alfven root in the inertial range has a small but non- 
vanishing 8B\\ which increases with k±, thereby causing 
the A m to decrease with k±. The VM solution for the 
Alfven root (black) A m arid dispersion relation are plot- 
ted against k±pi in Figures l2"al and [2bl for ft = 0.01, 0.1, 1, 
and 10 (dash-dotted, dotted, solid, and dashed respec- 



tively), ft = v^/v z A and 



Bq/ 'Airrriin is the Alfven 



speed. The solution assumes an inertial range spectral 
anisotropy as described by the Goldreich-Sridhar model, 
fen = feV 3 fe^ 3 , where ki is the isotropic outer-scale of 



the turbulence and is taken to be kipi = 10 , con- 
s istent with in sit u solar wind measurements at 1 AU 
Howes et aLll2008t) . 



Continuing along the critical balance cascade to kinetic 
scales naturally produces spectral anisotropy with k± ^ 
fen. So, we next consider the Alfven root with k±pi > 1 
and k\\di < 1. At these scales, the Alfven wave tran- 
sitions into the kinetic Alfven wave (KAW). The KAW 
is dispersive, damped, and much more compressible than 
the Alfven wave. The dispersive nature of the root steep- 
ens the magneti c energy spectrum to a = 7/3 in the 
unda mped case ( Howes et al.l |200§I : ISchekochihin et all 



l2009h . while t he inclusion of damp in g steepens the spec- 
trum further ( Howes et al.l l2011bl lcl; iTenBarge fc Howesl 



l2012a|) . Damping here refers to collisionless wave-particle 
interactions, primarily ion transit time damping on the 
pa rallel magneti c field that peak s at ion scales for ft > 
1 (|Barnesl Il966t iQuataertl fl998l) and electron Landau 
damping o n the parallel electric field that peaks at elec - 
tron scales (|Howes et al.ll2006tlSchekochihin et al.ll2009l) . 
The spectral anisotropy of the KAW is assumed to 



scale as fc || cx k { llljaitieil izuua iriowes et al.l l^l 
ISchekochihin et all 120091: ITenBarge fc Howedl20i2bi r 



1/3 



dGaltierl [20061: iHowes etHI 120081: 
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FIG. 2: (Color online) (a) The predicted variance anisotropy based on linear Vlasov-Maxwell theory for the Alfven 
(black), fast (red), and slow (blue) roots as a function of wavenumber. Green lines represent a fiducial solar wind 
mixture of the three modes composed of 90% Alfven, 9% slow, and 1% fast waves. Dot-dashed, dotted, solid, and 
dashed lines represent = 0.01, 0.1, 1, and 10 solutions respectively. The Alfven and slow roots are assumed to 

have a spectral anisotropy given by critical balance with an outer scale kipi = 10~ 4 and the fast root has Oub = 3° 

1/3 

up to kudi = 1 and fcii ex k j_ thereafter, (b) The dispersion relations of the same wave modes. 



The increased compressibility of the KAW can be seen 
at scales k±p t > 1 in Figure [2a] where the KAW A m 
is seen to have a f3i dependent plateau. The KAW A m 
in the dissipation range also depends upon the ion-to- 
electron temperature ratio, T,/T e . An analytical form for 
the KAW A m can be derived (see Appendix IA3[) within 
the frame work of electron reduced M HD (ERMHD) de- 
veloped in ISchekochihin et al. I (12009ft . 



2 + A(l+T e /T z ) 
Pi{l + T e /Ti) 



(1) 



The KAW A m has no wavenumber dependence and thus 
plateaus at a value determined by and Ti/T e . The 
^ and temperature ratio dependence of the KAW A m 
derived from VM theory is plotted in Figure El where 
the value for the A m is averaged across the plateau at 
k^ Pl e [3,4]. 

The behaviour of the Alfven root becomes more com- 
plicated for scales k^di > 1. When k^di > 1 and k±pi < 
1, the Alfven root becomes the Alfven ion cyclotron 
mode. This mode is strongly cyclotron damped and 
characterized by a left-handed magnetic hclicity and left- 
handed electric field polarization. Because of the strong 
damping and in situ solar wind observations suggest- 



ing this mode is less commo n than the KAW ( He et all 
l2012fc iPodesta fc Garvl[20rl . we do not consider it fur- 
ther. When k\\di > 1 and k±pi > 1, the behaviour of 
the Alfven root has not be en thoroug hly explored. The 
ion Bernstein wave (IBW) ([Stixl 119921 ) is conventionally 
assumed to couple to the KAW at the ion gyrofrcquency; 
however, the KAW root that exists for k^ dj < 1 may 
continue undamped at higher frequencies (|Howes et al.l 
120081: iKlein et al.ll2012at ISahraoui et al.ll201lh . Since the 
behaviour in this region of wavenumber space is uncer- 
tain, we do not attempt to describe it here. 



B. Slow Mode 

The majority of in situ solar wind observations of 
the inertial range suggest the component responsible for 
most of the measured free solar wind compressibility are 
pressure balanced structu res (PBSs) ([Bruno fc Carbona 
120051: iTu fc Marschlfl995h : however, linear PBSs are de- 
generat e with the fc || = 0, non-propagating limit of slow 
waves ([Kellogg fc Horburvl 120051: ITu fc Marschl Il995l) . 
For simplicity, we classify linear PBSs as slow modes. 
The association of the compressible portion of the solar 



4 




smaller than the Alfven cascade time, <C 1. 

For slow waves in the MHD limit, A m = k?Jkj_ (see 
Appendix I A II for the derivation of this equation), and the 
VM solution does not deviate significantly from the MHD 
solution until finite Larmor radius effects become domi- 
nant at k±pi ~ 0.1. The VM solution for the slow mode 

A m (blue) is plotted in Figure l2"al assuming a passive cas- 

2/3 

cade of slow waves with hug oc k± s — a more anisotropic 
cascade decreases the A m below that in the Figure but 
does not alter the qualitative behaviour. The A m for 
a PBS, i.e., fcii = slow mode, is identically zero since 
SB± = for these modes. 

Since the slow mode is strongly damped at large k\\ 
unless Ti <C T e , we will exclude a discussion of the be- 
haviour of this root for k\\dt > 1. 



FIG. 3: The ft dependence of the KAW root for 
Ti/T e = 10, 1, and 0.1 (dotted, solid, and dashed 
respectively) averaged across k±pi £ [3,4]. 



wind to PBSs is due to the measured anti-correlation of 
the t hermal and magnetic pressure ( Burlaga fc Qgilviel 
1 19T0h or an ti-correlation o f density and magnetic fi eld 
magnitude ( Roberts! Il990l : IVellante fc Lazarus _ 19871) at 
inertial range scales. Recent analyses ( Howes et al.1 
l2011aHKlein et al.ll2012bh exploring the more telling anti- 
correlation of density and parallel magnetic field suggest 
that the compressible portion of the solar wind is pri- 
marily composed of propagating slow modes rather than 
PBSs. These analyses suggest that on average Alfven 
modes comprise ~ 90% of the energy and slow modes 
~ 10% of the energy. 

The notion that the warm, Ti/T e ~ 1, ft ~ 1 
propagating slow mode exists in the solar wind defies 
conventional wisdom tha t the slow m ode is strongly 
damped in such plasmas ( Barnes! 19661). so we here elu- 



cidate this point. Theory ( Goldreich fe Sridharl 
ILithwick fc Goldreichl l200lj 



Schekochihin et al. 
Lazarianll2002l 



1997 



2009) 



2003) 



and numerical simulations (jCho 
suggest the slow mode does not have its own active 
turbulent cascade; rather, the slow mode is passively 
cascaded by the Alfvcnic turbulence. The slow mode 
damping rate is proportional to the parallel wavenum- 
ber, 7s cx k«g, and the strong damping of the slow mode 
suggests js/CJs — 1) where feiig, 7s, and ws are the par- 
allel wavenumber, linear damping rate, and frequency of 
the slow mode. However, the passive cascade of the slow 
mode implies that the slow modes are cascaded by the 
Alfven cascade on the Alfven timescale, u>a = k\\A v A- 
Therefore, at a given k^_, the parallel wavenumber of the 
slow mode compared to the parallel wavenumber of the 
Alfven mode determines the strength of the damping rel- 
ative to the cascade rate: at a particular scale, k*j_, a slow 
mode with k?< s <C k* A can be passively cascade before 



"\\A 



being damped since the slow mode damping rate will be 



C. Fast Mode 

At inertial range scales, the fast mode is well described 
by MHD, whose solution provides an A m identical to that 
of slow waves: A m = k?,/kj_. However, the distribution 
of fast waves in wavenumber space is not as well con- 
strained as that of slow modes since the fast mode is not 
strongly damped for parallel propagation. Further, com- 
pressible MHD turbulence simulations indicate that the 
fast mode is cascaded i sotropically in wavenumber space 
(ICho fc Lazarianl 120031) . Therefore, a turbulent cascade 
of fast modes has a A m that can take on all possible 
values, and the measured average value of the A m will 
depend sensitively on the distribution of fast modes in 
wavenumber space. To represent the fast mode turbu- 
lence, we plot in Figure [2a| the A m for a fast wave (red) 
with 6kB = arccos (fc ■ Bo/|fe||-Bo|) = 3°. We choose this 
value as representative because the average A m of an 
isotropic distribution of fast modes will be dominated 
by those modes with O^b — and the largest measured 
A m in a large en semble of solar wind data is ~ 500 
(ISmith et al.ll2006h . 

In the dissipation range at scales kndi > 1, the fast 
mode transitions to a parallel whistler (k±pi < 1) or 
an oblique whistler (fcj_Pi > 1) wave with w/Qj > 1. 
The whistler mode is wel l described by the electron 
MHD (EMHD) equations (jKingsep et al.lll990l) . which 
describes phenomena on scales kdi > 1. Note, one must 
take care when applying the EMHD equations, because 
EMHD is only valid for T < T e . For scales k±di > 1 
and kudi < 1, EMHD descri bes the cold ion limit of 
KAWs and not whistler waves (iHirose et al. ll2004l : lHowesl 
120091 : Eto et al. 1 12004 ISchekochihin et all I2009D . Solv- 
ing the EMHD equations provides A m = 1 + 2fcjj/fc^ 
(see Appendix IA 2[) , which is a good approximation of 
the VM solution and can again attain all values be- 
tween 1 and infinity and is sensitively dependent upon 
the wavenumber distribution. Therefore, if the propa- 
gation angle or whistler wavenumber distribution does 
not change from that of the inertial range fast modes, 



5 



the average A m will approximat ely double in the dissi 



pation range. However ; theory (ICho &: Lazarianl 



20041: 



iGaltier fc Bhattacharieeil2003HNarita fe Garvll201C ) and 



simulation (ICho fe Lazarianl |2004| iDastgeer et all 20001 : 
ISaito et "aT. 1 120101: ISvidzinski et al.ll2009h suggest the cas- 
cade of whistler waves is highly anisotropic in the same 

1 /3 

sense as critically balanced KAWs, kn oc fc ± . The 
anisotropic nature of the dissipation range cascade sug- 
gests the average A m of a distribution of whistler waves 
will decrease rapidly to a independent value A m ~ 1 
as the cascade progresses to smaller scales. The whistler 
portion of the fast modes plotted in Figure [2a] follows 
the critical balance prediction described above, with 

1 /3 

OkB = 3° up to k\\di = 1 and then fen oc kj_ . Note that 
the A m of whistler waves has a very weak fti and ion-to- 
electron temperature ratio dependence. The apparent Pi 
dependence of the whistler branch in Figure da] is due to 
the fast- whistler break point being at k\\di = k\\pi/ \/Jk\ 
therefore, the j3i dependence in the Figure would vanish 
if the x-axis were normalized to the ion inertial length. 

For completeness, the fast mode to whistler transition 
Am is shown for several different propagation angles, OkB, 
in Figure [4a] The A m agrees well with the predictions 
above, except in the O^b = 75° case. For such highly 
oblique angles, the fast mode transitions to an IBW: at 
scales k±pi > 1 and k\\di < 1, the fast mode transi- 
tions to the electrostatic IBW with frequency approxi- 
mately equ al to an integer multiple of the ion cyclotro n 
frequency (|Howesl l2009t iLi fc Habball l200lt Istixl Il992t) . 
The transition to the IBW is clear in the fast/whistlcr 
wave dispersion relations in Figure I4bl corresponding for 
the same angles as the A m presented in 0a) Although the 
IBW is electromagnetic during the transition from elec- 
tromagnetic f ast or Alfv en modes, IBWs are dominantly 
electrostatic ( Stixl fl992l) . The properties of the transi- 
tion electromagnetic IBW have not been fully explored 
in the literature, so we will not consider the potential 
contribution of IBWs to the dissipation range A m . 

The wavenumber distribution of whistler waves gen- 
erated in a local instability such as the parallel firehose 
( Gary et al.lll998h and whistler anisotropy instabilities is 
not easily determined and depends upon the instability, 
so we will not consider them further. 



III. COMBINATIONS OF MODES AND A m 

Assuming the solar wind consists of combinations of 
the three wave modes described in the preceding section 
both in the inertial and dissipation ranges, we need to 
consider the effect of superposing the modes. The three 
linear modes can be combined into a total A m measure 



A m T = 

AC C ±A + (1 



AC)[FS C ±F + (1-FS)C ±S ] (2) 



where AC is the fraction of Alfven to total energy, FS is 
the fraction of fast to total compressible (fast plus slow) 
energy, C± = (\5Bj_\ 2 ) / (\5B\ 2 ) , and subscripts A, F, and 
S refer to Alfven, fast, and slow modes. Note that AC £ 
[0, 1] and FS € [0, 1]. C\\ and C± for each mode can both 
be written in terms of the A m as 



Cn = 



f 



and 



1 + A n 



•Am 
l + Ar, 



(3) 



(4) 



A. Inertial Range 

The assumptions that only slow modes with fc|| 3> 
k± are weakly damped and approximately parallel fast 
modes dominate the inertial range fast mode A m lead to 
asymptotically small and large values of the slow and fast 
mode A m s, respectively. Due to the asymptotically large 
and small values of the component A m s, we can estimate 
the value of A m T in the inertial range provided AC =/= 1 
and FS ^ 1 by using the A m of each mode from Figure 
I2al to estimate the compressibilities, which have only a 
very weak j5i dependence: C±a — 1, C\\a <C 1, C±f — 1, 
C\\ F <C 1, C±s <C 1, and Cng ~ 1. Therefore, in the 
inertial range, 



A 



mT 



AC + (1-AC)FS 
(1 - AC)(1-FS) ' 



(5) 



AC C\\ A + (1 - AC) [FS C\\ F + (1 - FS)C\\ S ] 



Note that whether the slow modes are in fact propagating 
slow modes or PBSs does not effect the estimate because 
C±s is asymptotically small in either case. Therefore, the 
inertial range Am depends only on the ratios of Alfven 
to total energy and the fast to total compressible energy 
ratio and is unable to differentiate between propagating 
slow modes and PBSs. Note that for observe d average so- 
lar wind values of AC ~ .9 and FS < 0.1 ( Howes et al.l 
l2011al: iKlein et al.l l2012bl) . the A m can be further re- 
duced to A rn T — AC I (1 — AC). We take as fiducial solar 
wind values AC = 0.9 and FS = 0.1, or cquivalently 90% 
Alfven, 9% slow, and 1% fast wave energy. Although we 
take FS = 0.1, there is very little quantitative difference 
between FS = 0.1 and (see Figured]). 

Plotted in green in Figure [2a] is a single mixture of 
the three modes representing fiducial solar wind values of 
AC = 0.9 and FS = 0.1. Clearly, the A m has virtually 
no Pi or wavenumber dependence in the inertial range, 
as expected from the above analysis. To explore the de- 
pendence of the A m on fluctuation composition, we plot 
in Figure [5] the A m T given by equation ([5]). To confirm 
the quality of the estimate given by equation ([5]) , we plot 
in Figure [6] the A m T given by equation ([2]) (thick) and 
the estimate given by equation ([S]) (thin). The Figure 
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FIG. 4: (Color online) (a) The predicted variance anisotropy based on linear Vlasov-Maxwell theory with = 1 and 
Ti = T e for the fast/whistler root for different propagation angles, QkB = 3°, 15°, 30°, 45°, 60°, 75° (long-short 
dashed, long dash-dotted, short dash-dotted, long dashed, short dashed, dotted), and 3° transitioning to critically 
balanced whistler (solid). When the QkB = 75° root transitions to an IBW, the color changes to red. The 8kB = 3° 
root terminates at k±pi ~ 3 due to the solver's inability to resolve the high frequency of the root, (b) The dispersion 

relation for each root. 



demonstrates that equation ([5]) provides a good estimate 
for the total A m in the inertial range. 

Also highlighted by the Figure is the degeneracy of 
the A m , since different mixtures of modes can replicate 
nearly identical A m s in the inertial range. This implies 
that the A m when used alone is not a good metric for 
differentiating modes in the inertial range; however, the 
Am can be a useful secondary metric. For instance, when 
used together with density-parallel magnetic field corre- 
lations to identify the proportion of fast to slow wave 
energy, the A m can provide an estimate of the Alfven to 
compressible energy proportion, thereby quantifying the 
total population of the solar wind. 



B. Transition and Dissipation Range 

The behaviour of the A m in the transition between 
the inertial and dissipation ranges and in the dissipation 
range is markedly different from the inertial range. Un- 
like the inertial range, where the asymptotically large and 
small values of the Alfven, fast, and slow waves leads to 
an A m that is controlled by the mode fractions AC and 
FS, the KAW dominates the dissipation range A m for 



typical solar wind parameters (see Figure l2al for fcj_P« > 1, 
where the A m for the mixture closely follows the A m for 
KAWs). 

The strong /3, and Ti/T e dependence of the KAW 
A m implies that for certain values, namely > 1, the 
whistler and KAW dissipation range A m become approx- 
imately degenerate (see Figure [2a| for k±pi > 1). How- 
ever, the behaviour through the transition range for the 
Alfven to KAW and fast to whistler mode transitions 
differs considerably. For example, the curves in Figure 
ED with (AC,FS) = (0.9,0.1) (solid black) and (0.1,0.9) 
(short dash-dotted black) are nearly identical in the iner- 
tial and dissipation ranges, but exhibit very different be- 
haviour across the transition range. For this reason, the 
A m is best presented as a function of wavenumber rather 
than as a quantity averaged over a band of wavenumbers, 
as has often been done i n the literature ( Hamilton et al.l 
mith et al. l l2006ll2012D . 



IV. MEASURING A m 

The conventional method for calculating the A m is as 
defined in $2 however, different methods of averaging can 
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FIG. 5: Contours of the inertial range log 10 (A m T) for 
different mixtures of fluctuations determined by the 
fractions of Alfven to total energy (AC) and fast to 
total compressible energy (FS). 



be performed. We here consider two physically motivated 
alternatives for measuring the A m and compare them 
to the linear prediction. We also explore the validity of 
using linear theory to describe the A m in a fully nonlinear 
turbulent situation. 




FIG. 6: The sum of the Vlasov-Maxwell Alfven, fast, 
and slow modes for = 1 (thick) and the inertial range 

estimate given by equation (thin) for different 
combinations of the three roots. The different line types 
represent Alfven to total energy (AC) and fast to total 
compressible energy (FS) fractions: 
(AC,FS) = (0.9,0.1), (0.99,0.1), (0.8,0.1), (0.6,0.1), 
(0.1,0.9), and (0.5,0.99) are solid, dotted, short dashed, 
long dashed, short dash-dotted, and long dash-dotted. 



A. Synthetic Solar Wind Measurements 



In the solar wind, the mean magnetic field direction 
is constantly changing and can sweep through a wide 
range of angles over a period of tens of minutes. As 
such, defining 8B\\ = Bq ■ SB becomes questionable since 
Bo is averaged over long, global, periods relative to the 
rapidly fluctuating field — this is especially problematic 
when measuring in the dissipation range when 5B\\ is 
even smaller and more rapidly fluctuating. The poor def- 
inition of parallel and perpendicular in the global anal- 
ysis will pollute the parallel energy with perpendicular 
energy, thereby decreasing the measured solar wind A m - 
Therefore, a local analysis employ i ng a local mean mag- 
netic field (ICho fc yishniad l2000i: jHorburv et all 120081: 



M aron fc Goldreichl 1200 ll iPodestal I2009T) must be em- 
ployed to accurately measure cVBm . The local analysis 
has the added advantage of being able to differentiate 
between different 6vb measurements, where 9vb is the 
angle between the mean magnetic field and the solar wind 
flow velocity. This can be helpful because the A m will 
have different incrtial/dissipation range breakpoints for 
KAWs and whistlers when plotted against fcy or k±. 

The conventional A m is constructed by calculating sep- 
arately the perpendicular and parallel fluctuating mag- 
netic energies and finding their quotient, 



A — 



(\SB A 



(\SB\\\ 



(G) 



However, a similar measure could be constructed from 
the normalized perpendicular and parallel energies, 
(\SB ± \y\SB\ 2 ) and (\6Btf/\5B\*), 



A — 



(\SB ± \y\6B\ 
(\5B l{ \y\5B\- 



(7) 



This measure has the physically motivated advantage of 
avoiding possible small denominators due to small SBu. 
Another possible measure of the A m is 



A m = (ISB^/ldB. 



(8) 



This method has the conceptual advantage of maximiz- 
ing any local effect of the mean magnetic field since 
it averages the A m at each point rather than sepa- 
rately averaging the energies. However, the expression 
(|(5-Bj_| 2 /|(5-B|| | 2 ) will be dominated by those terms with 
small denominators, which could make it an unphysical 
measure of the A m 

To explore the effect different averaging procedures 
have on the A m , we employ synthetic spacecraft data. 
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FIG. 7: Synthetic spacecraft data populated with 90% 
Alfven, 9% slow, and 1% fast wave energy averaged over 
k±p.i € [0.003,0.03]. The three different lines represent 
different averaging procedures for calculating the A m - 



The details for c onstructing t he synt hetic spacecraft data 
are discussed in iKlein et al.1 (l2012bh . so we here only de- 
tail the relevant parameters. A three-dimensional box 
is populated with a spectrum of all three wave modes, 
where the proportion of each mode is given by AC = 0.9 
and FS = 0.1, corresponding to 90% Alfven, 9% slow, 
and 1% fast wave energy. The Alfven and slow modes sat- 
isfy critical balance, with all modes less than the critical 
balance envelope equally populated. The fast modes are 
isotropically populated from 8kB — 5 — 85°, with equal 
energy at each angle. To represent solar wind turbulence, 
the phase of the wave modes is randomized. The data 
is sampled by adverting the turbulence past a stationary 
"spacecraft" with velocity v = lOw^ at a fixed angle with 
respect to the mean magnetic field and fixed sampling 
rate. 

Figure [7] presents the three different methods for cal- 
culating the A, n . The A m is averaged across the interval 
k ±Pl e [0.003,0.03] with 6 vB = 85°, where 9 vB is the 
sampling angle analogous to that in solar wind. Sam- 
pling angle, 9 v b, was varied from 5 — 85° and found to 
have no significant effect on the A m in the inertial or 
dissipation range, but sampling angle does alter the be- 
haviour in the transition region. 

For AC = 0.9 and FS = 0.1, the predicted A m from 
linear theory is A m T — 10. From the Figure, it is clear 
that the conventional definition of the A m , equation (|6|), 
(solid) computed from a spectrum of wave modes agrees 
best with linear theory and is the one we will continue to 
employ. The averaged anisotropy, equation (JT]) , (dotted) 
approach is a poor measure of the A m because at any 
given point in space, the superposition of a collection 
of wave modes can lead to anomalously small values of 
5B» due to cancellation. While the definition based upon 



normalized energies, equation ([8]), (dashed) differs from 
the linear prediction, it captures the correct qualitative 
behaviour and may be a safer measure to use in some 
circumstances. 



B. Nonlinear Simulation Measurements 

Here we present results from a collection of fully non- 
linear gyrokinctic turbulence simulations performed with 
/3i = 0.1 and 1, Ti/T e = 1, and a realistic mass ra- 
tio, mi/m e — 1836. The simulations were performed 
using the Astrophys ical Gyrokinetics Code, AstroGK 
( Nu mata et al ] I2010D . All runs presented herein are 
driven with an oscillating Langevin antenna coupled to 
the p arallel vector potentia l at the simulation domain 
scale (jTenBarge et al.l |2012() . Relevant parameters for 
the simulations arc given in table HI where e = Pi/ Lq <C 1 
is the gyrokinetic expansion parameter, Aq is the an- 
tenna amplitude, and v s is the collision frequency of 
species s. The expansion parameter sets the parallel 
simulation domain elongation and is determined by as- 
suming critical balance with an outer-scale kipi = 10~ 4 : 



* 5 /3 

k/ )/(l + k± ), where k = kpi and sub- 



,-1/3,-2/3,-, 

script naught indicates simulation domain scale quanti- 
ties. The antenna amplitude is chosen to satisfy critical 
balance at the domain scale so that the simulations all 
represent critically balanced, strong turbulence. 

An example instantaneous one-dimensional perpen- 
dicular magnetic energy spectrum composed by over- 
laying the four f3i = 1 simulations is plotted in 
Figure [5] — the spectrum is produced via conventional 
Fourier analysis techniques. The inertial range simula- 
tion (red) has a spectral index ~ —3/2, which steep- 
ens (green) to —2.8 (blue) before rolling off exponen- 
tially (cyan) as the electron collisionless dissipation be- 
comes increasingly strong. The dissipation range sim- 
ulations /31D and / 31ED have been explored in detail 
(jHowes et all201ld : lTenBarge fe Howeil2012af) . The full 
spectrum agrees well with recent solar wind observations 
(jAlexandrova et al.ll201l[ ). 

Figures l9al and [9bl overlay the A m for the /3, = 0.1 and 
1 AstroGK simulations. The A m is produced by comput- 
ing, via conventional Fourier techniques, the perpendic- 
ular and parallel energy with respect to the global mean 
magnetic field. D ue to the restrictions of gyrokinetics 
( Howes et al.ll2006l ) and the method of driving, the simu- 
lations are almost purely populated by Alfven / KAWs. 
Therefore, we plot as a dashed line in each Figure the 
VM linear solution for the Alfven root. The agreement 
between the linear prediction and the nonlinear simula- 
tion is excellent up to the point of strong damping at 
k±Pi ~ 10, where agreement with linear theory is ex- 
pected to breakdown. The good agreement supports the 
validity of using linear theory to describe the A m of fully 
nonlinear turbulence. 

Note that in constructing Figures [9a| and [9b] a global 
mean magnetic field has been used. This is justified, 
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TABLE I: Parameters for AstroGK simulations used in Figures [51 [5al and [9b] 
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FIG. 8: (Color online) Composite one-dimensional 
perpendicular magnetic energy spectrum composed of 

four overlapping /3j = 1 AstroGK turbulence 
simulations spanning from the inertial range to deep 
into the dissipation range. 



because in typical numerical simulations, the global mean 
field direction is fixed and \5B\/B -C 1. Although a 
local mean field can be defined in the same sense as in 
the solar wind, the global mean field provides a relatively 
accurate measure of the A m in numerical simulations of 
turbulence with a fixed global mean field direction and 
\SB\/B « 1. 



V. SOLAR WIND MEASUREMENTS 



We here discuss some of the recent solar wind anal- 
yses of the Am, present new solar wind measurements 
of the Am, and compare the new measurements to our 
predictions. 



ISmith et al.l ( 20061 ) presents the most comprehensive 
collection of the solar wind inertial range A m made to 
date, incorporating 960 data intervals recorded by the 
ACE spacecraft. The data is divide d into periods of ope n 
field lines and magnetic clouds ( Burlaga et al.l Il98l[) . 
They find that the A m is proportional to a power of 
Pi or \8B\/Bq but are unable to identify which is the 
source of the functional relationship since Pi and \SB\ 
are themselves posit ively correlated in the solar wind 
( Grappin et al.lll990l ). They follow the identification of 
the A m relationship with a discussion of possible sources 
for the relationship based on a variety of considerations, 
some of which con t radict the discussion in fJTT] of this 
paper. I Smith et~aT. I ([20061) do not consider the superpo- 
sition of modes that has been shown to exist in the solar 
wind, neglect the contribution of slow modes due to the 
belief that they are completely damped outside of local 
excitation, consider Alfven waves to be completely in- 
compressible (A m — oo), consider PBSs but assert that 
the j3i dependence of 6B\\ of PBSs could contribute the 
the /3i dependent Am despite A m PBS = for all Pi. 

The findings of ISmith et al. (2006) seem to be in di- 
rect contradiction to the conclusions drawn in §111 A[ 
however, there are possi ble reasons f o r the apparent dis- 
crepancy. As noted in ISmith et al.l ( 2006[ ) , the proton 
temperature and fluctuating magnetic field are positively 
correlated. It has also been shown that the proton tem- 
perature and bulk solar wi nd speed are positively cor - 
related (|Burlaga fc Qgilvid Il973t iNewburv etall Il998h . 



Therefore, the measured A m might depend on any or 
a combination of Pi, \5B\, or Vsw- If the underlying 
dependence is actually on vsw, then the observed vari- 
ation of the A m could be from different types of solar 
wind launched from different regions of the sun having 
a variable population of Alfven to compressible compo- 
nents. For instance, slower wind may have evolved to a 
more Alfvenic state, because fast modes would have had 
time to be dissipated in shocks and slow modes could 
have collisionlcssly damped, leaving a primarily Alfvenic 
and PBS population. Since no meas ure other than A m 
binned by Pi or \5B\/Bq was used in ( Smith et al.ll2006T ). 
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FIG. 9: (Color online) A comparison of the linear Vlasov-Maxwell solution (dashed) to overlaid nonlinear 
gyrokinctic simulations of Alfvcnic turbulence spanning from the incrtial range to deep into the dissipation range for 

(a) ft = 0.1 and (b) ft = 1. 



it is difficult to draw a firm conclusion. 

lHamilton et al.1 ( 2008T) extend the s tudy of the same 
dataset employed in lSmith et al.l ( 20061 ) to the dissipation 
range, defined to be 0.3 Hz < / < 0.8 Hz . They find 
that dissipation range fluctuations are less anisotropic 
th an those in the inertial range. Although not discussed 
in lHamilton et al.l ( 20081) , this result is consistent with 
the mixtures of wave modes discussed in this paper. The 
linear predictions of KAW A m presented in Figure [3] for 
Ti/T e = 1 pass through the core of the measured A m 
for open magnetic field data for the measured A m in the 
dissipati on range presented in the lower panel of Fig- 
ure 8 of lHamilton et al.l (|2008l ). Further, the KAW A m 
for TjT e = 0.1 and Ti/T e = 1 bound the core of the 
low ft magne tic cloud data fo r the m easured dissipation 
range A m in lHamilton et all ( 20081 ). T he temperature 
ratio is not provided for the dataset in lHamilton et al.l 
(120081) : however, magnetic clouds are typically character- 
ized by low proton tem peratures ( Osherovich et alJll993t 
iRichardson et al. 1997), so T j /T K < 1 for the cloud data 
is expected. iGarv fc Smithl ( 20091 ) also noted that the 
ft > 0.1 portion of the same dataset is well fit by the 
KAW solution; however, the temperatu r e rati o depen- 
dence was not explored. IGarv fc Smithl ( 20091 ) suggest 



the shallow slope at low ft might be suggestive of a 
whistler component. Another possible explanation for 
some of the very low v alues of the dissipation range A m in 
lHamilton et al 1 (I2008D is the use of a global measurement 



of the magnetic field, which will tend to decrease the mea- 
sured A m , especially in the dissipation range. Also, due 
to the limited high frequency information available from 



ACE measurements fnyquist — 2 Hz (S mith et al.lll998l) . 
the region defined as the dissipation range corresponds to 
the transition region between the inertial and dissipation 
ranges. 

ISmith et al.l ( 20121 ) return again to the same ACE 
dataset but focus on ft > 1 intervals and augment the set 
with 29 additional high ft ACE measurements. As such, 
this dataset suffers from the same high frequen cy limita- 
tions noted above f or the lHamilton et all ( 20081 ) analysis. 
ISmith et al.1 ( 20121 ) suggest that the measured dissipa- 
tion range A m is consistent with KAWs, but the iner- 
tial range A m cannot be explained by a purely Alfvenic 
population. This observation is partially used to con- 
clude that the Alfven/KAW model cannot explain solar 
wind measurements. However, their A m measurements 
agree well with a solar wind population dominated by 
Alfvcn/KAWs with a small component of slow waves. 
Therefore, we find no contradiction between the ACE 
A m measurements and the Alfven mode dominant A m 
model constructed herein. 

Using STEREO measurements, iHe et al.l (|2012t ) em- 
ploy a variant of the A m , SBu /SBj_ , as a secondary metric 
to the magnetic hclicity to differentiate between KAWs 
and whistlers, which both have right-handed hclicity for 
oblique propagation (|Howes fc Quataertl[2"010f ). However, 
they consider fast modes propagating at fixed 9kB = 60°, 
80°, and 89°, and only the Bub = 60° root connects to 
an oblique whistler. The other two oblique roots connect 
to IBWs, which is made cl ear in the upper right panel of 
Figure 3 of lHe et all ( 20121 ). where the dispersion relation 
for these two roots do not extend above uj/Cli = 1. Due 
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to this mischaracterization of the roots, they mistakenly 
find the oblique whistler root A m < 1. As noted in mi CI 
the whistler (fast mode for k\\di > 1) A m > 1 for all prop- 
agation angles and Pi and is thus of the same orientation 
as KAW and does not provide a useful secondary metric 
to the helicity. 

Using Cluster measurements, ISalem et al.1 (120121) use 
magnetic compressibility, which is related to the A m in 
equation ([3]), as a secondary metric to the ratio of the 
electric to magnetic field, |<5-E|/|<5.B|, which is degener- 
ate between KAWs and whistlers. Although the dissipa- 
tion range asymptotic behaviour of the fixed propagation 
angle KAWs and the whistlers at certain angles is simi- 
lar, the transition from the inertial range to dissipation 
range for the two modes is different, and it is the tran- 
sition region that i s used to differentiate the two modes. 
ISalem et al. (2012) find the Cluster data to be most con- 
sistent with a spectrum of KAWs. 

Although much of the work does not directly address 
the A m , some recent solar wind analyses of t he magnetic 
field have moved toward three dimensions (IChen et al.1 



l2011at iPodesta fc TenBared l2012t IWicks et all I2012D 



Presenting the full three dimensional structure of the 
magnetic field is of obvious potential benefit; however, 
there are complications due to geometrical and sampling 
effects of the fluctuations being advected pas t the space- 
craft ([Turner et al.l 1201 it I Wicks et al.l 1201 2h. To eluci- 
date the complications, wc follow iTurner et al. I (I20TI 
and choose a coordinate system with Bq in the z direction 
and V sw m the y direction. Under Taylor's hypothesis, 
the measured components of magnetic energy in the per- 
pendicular plane are then 



E x ,y(cj) = — / d 3 k\SB x J 2 6{co - k ■ V sw ). 
If a scaling between fcii and kj_ is assumed, then 



E x {lo) = J d 2 k ± E 2 D(k±)S(u; 



kyVsw ) cos (</>) 



E y (w)= / d 2 k ± E 2D (k ± )6(u; - k y V sw ) sin 2 {4>), 



where E 2 u(k±) = Ck ± ^ 1 is the two-dimensional spec- 
trum, £ is the one-dimensional sca ling exponent, and 
tan ((/)) = k y /k x . ITurner et all (|201lh find that E y /E x = 
l/£, implying a nonaxisymmetric energy distribution in 
the perpendicular plane. This complication is obviated 
by measuring A m in the standard way, since 

E±{u) = E x (lu) + E y (uj) = 



d 2 k±E 2D (k 1 _)8(uj - kyV s , 



avoids the angular sampling issue altogether. 



B. New Measurements 

To compare the analysis method outlined in ^ITI and HTT1 
to solar wind measurements, we choose a 5 day interval, 
2008 Feb 12 06:00 to Feb 17 06:00, of magnetic field data 
captured by the Stereo A spacecraf t. This interval is one 
of 20 similar intervals analysed by IPodesta fc TenBared 
(j2012h . where details concerning the data selection and 
analysis can be found. So, we here summarize only the 
most pertinent aspects of the measurement and analy- 
sis. The interval was chosen because it represents an 
interval of high-speed wind, which typically satisfies the 
assumption that the solar wind flow direction is in the 
heliocentric radial direction, R. Also, the interval is suf- 
ficiently long to include a statistically large sample of pe- 
riods during which the mean magnetic field, Bo, is nearly 
perpendicular to the radial direction, which in practice 
includes the range 84° < 9vb < 96°. Focusing on these 
orthogonal periods is helpful because it facilitates easier 
comparison to theory since the measured wavenumber in 
the solar wind corresponds more closely to the perpen- 
dicular wavenumber used in linear theory. 

The magnetic field data was analysed using wavelet 
techniques to determine t he local mean magnetic field, 
Bq(t ), at a given scale r ( Horburv et al.ll2008l : IPodestal 
I2009D . The parallel and perpendicular magnetic en- 
ergies are then given by SB^(t) 2 = 
and 5B ± (t) 2 = \5B(t)\ 2 - SB^t) 2 
B (t)/B (t). 

The proton plasma parameters for the interval were 
supp lied by the PLASTIC instrument ( Galvin et al.l 
120081 ). The relevant proton plasma parameters are V p = 
657 km/s, n p = 2.6 cm" 3 , T p = 1.7 x 10 5 K, (3 p = 0.65, 
where all quantities are averaged over the entire in- 
terval except j3 p whose median is taken because of its 
high variability. Electron thermal data is unavailable 
due to problems with the IM PACT Solar Wind E lec- 
tron Analyzer aboard Stereo ( Fedorov et al.l l201ll ) , so 
we assume the electron temper ature to be the averag e 
for fast wind streams found by INewburv" et all (11998), 
T e = 1.4 x 10 5 K. This assumption for the electrons 
implies T p /T e ~ 1 for the data interval. 

To convert spacecraft-frame frequency, f sc , to 
wavenumber, we employ Taylor's hypothesis, which 



B (t) -5B(t)\ 2 
where Bq(t) = 



states that 2nf s 



OJ, 



k ■ V sw , where 



the plasma rest-frame frequency. Fo r oblique Alfvenic 
and KAW fluctua tions, uj p <C fc ■ V sw ( Howes et al.ll20l3 : 
ISmith et al.ll2012l ). and we can assume 27r/ sc = kV sw . To 
compare to linear theory, we convert spacecraft frame fre- 
quency to wavenumber using kpi = 2ir Pif S c/V sw , where 
Pi is computed from the average solar wind proton quan- 
tities for the measured interval. 

In Figure [TUl we compare the measured A m for this 
interval (blue) to the A m from linear theory with /3.; = 
0.65 and Ti/T e = 1 for two different mixtures of wave 
modes that reproduce the solar wind A m in the inertial 
range: a mixture that is dominantly Alfvenic (black) with 
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(AC,FS) = (0.9,0.1) and a mixture that is dominantly 
fast mode (red dashed) with (AC, FS) = (0.1,0.9). Al- 
though both mixtures of linear wave modes reproduce 
well the behaviour in the inertial range, the fast mode 
dominant construction displays markedly different transi- 
tion and dissipation range behaviour; whereas, the Alfven 
dominant mixture fits well the slope of the transition 
range and the asymptotic value in the dissipation range. 
From this, we can conclude that the measured interval 
is most likely dominated by Alfven wave fluctuations at 
k±pi < 1 and KAW fluctuations for k±pi > 1. 

We speculate that the positive slope in the measured 
inertial range A m may be attributable to a decreasing 
fraction of slow wave energy due to collisionless damp- 
ing as the turbulent cascade proceeds to kinetic scales. 
Also, the approximate factor of two difference between 
the predicted A m from linear theory and the solar wind 
values in the transition and dissipation ranges arc likely 
due to averaging ft and pi over the entire 5 day inter- 
val. The A m of KAWs has a moderately strong ft de- 
pendence (see Figures l2"al and 15)). For the purposes of 
comparison, we use a median ft = 0.65. However, the 
A m of KAWs increases with decreasing ft, so periods of 
low ft will tend to dominate the average A m and shift 
it upward from the prediction based on a median value 
of ft = 0.65. Similarly, we use a value for pi based on 
quantities averaged over 5 days of data to determine the 
relationship between frequency and wavenumber. The 
variance associated with the averaged pi could cause an 
apparent horizontal shift of the data in the transition and 
dissipation ranges. 

For comparison, we also plot in Figure [10] the mea- 
sured Am for the Stereo A data computed via a more 
conventional global analysis (green). The 5 day inter- 
val of data was segmented into 120 1 hr intervals. In 
each 1 hr segment, the mean magnetic field was calcu- 
lated, the data rotated into mean field coordinates, and 
a Welch windowed FFT was applied to compute the par- 
allel and perpendicular magnetic energies. The parallel 
and perpendicular magnetic energies were then averaged 
over the 120 1 hr segments and binned into logarithmic 
wavenumber bins to smooth the spectra. The global A m 
(green) was then calculated from the binned and aver- 
aged magnetic field data. As expected, this method is 
inappropriate to compare to linear theory and produces 
a significantly reduced A m compared to the wavelet anal- 
ysis because the parallel magnetic energy is polluted with 
perpendicular energy, as discussed in §IV Al 

This analysis of solar wind data highlights the impor- 
tance of performing a local analysis and the weaknesses 
and strengths of the A m for determining the solar wind 
fluctuation composition: The A m is a poor tool to use in 
the inertial range due to degeneracies of different wave 
mode mixtures. However, if the inertial range A m is aug- 
mented with a second measure that identifies the the fast 
to total compressible energy fraction (FS), such as the 
density-parallel magnetic field correlation, the A m can 
provide a measure of the Alfven to total energy fraction 
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FIG. 10: (Color online) A comparison between the 
measured Stereo A A m data computed via a local 
wavelet analysis (blue) in the 84° < 9vb < 96° bin, the 

measured Stereo A A m computed with a global, 1 hr 
averaged, mean magnetic field (green), and the summed 
Vlasov-Maxwell A m with ft = 0.65 and Tj/T e = 1 
employing (AC,FS) = (0.9,0.1) (black) and 
(AC,FS) = (0.1,0.9) (red dashed). 



(AC). Also, the slope of the A m transition range and 
the asymptotic value in the dissipation range arc useful 
for identifying the solar wind fluctuation composition. 



VI. CONCLUSIONS 

We have developed a framework for interpreting the 
measured solar wind magnetic variance anisotropy (Am) 
and comparing the measurements to linear theory. In 
CHI we reviewed the linear properties of the collision- 
less counterparts to the three MHD wave modes, which 
are expected to be the primary constituents of the so- 
lar wind in those relevant regions of wavenumber space 
spanning the inertial and dissipation ranges where colli- 
sionless damping rate is small compared to the rate non- 
linear energy transfer. The A m of each wave mode for a 
range of plasma betas, ft, is shown in Figure l2al versus 
perpendicular wavenumber. 

In §1111 we examined how superpositions of the three 
wave modes affects the A m - Due to the asymptotically 
large A m of the Alfven and fast modes and asymptoti- 
cally small Am of the slow mode in the inertial range, 
the inertial range A m is dictated by the fraction of each 
mode rather than the individual behaviour of each mode 
and has little ft dependence. Note that linear pressure 
balanced structures (PBSs) are equivalent to fcy = 0, non- 
propagating slow modes, and are thus classified as slow 
modes. An estimate of the inertial range A m is given by 
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equation ([5]) , where the A m is determined by the fraction 
of Alfven to total (Alfven plus fast plus slow) energy and 
the fraction of fast to total compressible (fast plus slow) 
energy. For /?.; > 1, the dissipation range value of the 
A m for kinetic Alfven waves (KAWs) and whistler waves 
is approximately degenerate; however, the behaviour of 
the A m through the transition to the dissipation range 
differs in form for KAWs and whistlers. The degeneracy 
in the dissipation range highlights the value of displaying 
the A m as a function of wavenumber. 

In ^IV Bl we employ a suite of fully nonlinear gyroki- 
netic turbulence simulations that span from the inertial 
range to deep into the dissipation range for = 0.1 and 
1 to demonstrate that the predictions of the Am from lin- 
ear theory are applicable to nonlinear turbulence. The 
comparison between the linear prediction for the A m of 
Alfven waves and the A m measured in the nonlinear tur- 
bulence simulation are presented in Figure 9, where excel- 
lent agreement is seen up to kinetic scales where electron 
dissipation becomes significant, k±pt ~ 10. 

Some of the recent solar wind measurements of the A m 
and their interpretations are presented in W Al Since the 
A m of previous solar wind studies typically bins the A m 
over a band of wavenumbers and computes the A m via a 
global analysis, the data is difficult to interpret and com- 
pare to linear theory; however, the Am observed in previ- 
ous studies is mostly consistent with solar wind composed 
primarily of Alfvenic fluctuations and conforms well to 
the theory outlined in this paper. 

WBI presented a measurement of the solar wind A m 
from Stereo A performed using the methods outlined in 
< HV Al computed locally via wavelet analysis techniques, 
and plotted as a function of perpendicular wavenum- 
ber. We find excellent agreement between linear the- 
ory and the local solar wind measurement, which sug- 
gests that this interval of solar wind data is composed 
of approximately 90% Alfven, 10% slow, and negligible 
fast wave energy across the sampled range of scales span- 
ning the inertial range to the beginning of the dissipation 
range. Here, slow wave energy encompasses both PBSs 
and propagating slow modes because the A m alone can- 
not differentiate between the two. This analysis high- 
lights the advantages of computing the A m locally with 
a wavelet analysis and plotting the A m as a function 
of wavenumber rather than binning across wavenumber 
bands because the transition range breaks the degener- 
acy of the inertial range. Alternatively, a second measure 
such as the density-parallel magnetic field correlation can 
be used to break the degeneracy and ascertain a value for 
fast to total compressible energy fraction (FS) . The A m 
can then be used to measure the Alfven to total energy 
energy fraction (AC). 

In summary, we have outlined the salient features of 
linear theory as they relate to the magnetic variance 
anisotropy, provided a procedure for producing and com- 
paring solar wind measurements to predictions from lin- 
ear theory, demonstrated the validity of this technique 
via both nonlinear turbulence simulations and solar wind 



measurements, and found previous studies of the solar 
wind magnetic variance anisotropy and new solar wind 
data to be consistent with a dominantly Alfvenic popu- 
lation in the inertial and dissipation ranges. 
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Appendix A: Fluid Limits of the Magnetic Variance 
Anisotropy 

Here we provide analytical derivations for the A m 
from fluid theories. Without loss of generality, we will 
assume an equilibrium magnetic field JBo = BqZ and 
k = (k±, 0, fc||) for all of the systems considered. For the 
purposes of constructing the A m , we need only consider 
the eigenfunctions of the magnetic field. 



1. Magnetohydrodynamics 

The Magnetohydrodynamic (MHD) equations can be 
written as 



| + v.(,«) = o, 



(Al) 



dt \p\ 



cE + u x B = 0, 
V • B = 0, 

dB 

— + cV x E = 0. 
at 

The linear dispersion relation for this system is 

(w 2 - v 2 A k 2 ) L 4 - oj 2 k 2 {c 2 + v 2 A ) + k\k 2 v\c 2 s 
0, 



(A3) 

(A4) 
(A5) 

(A6) 
(A7) 
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where the first term corresponds to the Alfven root and 
the two remaining roots correspond to the fast and slow 
compressible roots. Finally, the eigenfunctions for the 
system can be expressed as 



—k\h 



- v A k \\)-n~ = 

u 2 (c 2 + v 2 A )-c 2 s v 2 A k 2 SB 



(A8) 



and 



5B n 



Bo 



(A9) 



The Alfven root corresponds to uj 2 = v A k 2 , which im- 
plies &B|| = 0. Therefore, the A m for the MHD Alfven 
root is formally infinite. 

The fast and slow compressible roots correspond to 



7^ v A k 2 , which implies SB y = 0. Therefore, 8B\ 
5B 2 and after manipulation equation (|A8j) reduces to 



3. Electron Reduced MHD 

Electron Reduced MHD ( ERMHD) was introduced 
in ISchekochihin et all ( 20091 ) and is the k±pi >• 1, 
m e /mi <C 1 limit of gyrokinctics. The system general- 
izes the EMHD equations for low-frequency, anisotropic 
(fc|| -C k±) fluctuations without assuming incompress- 
ibility, and kinetic Alfven waves are described well by 
ERMHD. The equations of ERMHD are most simply ex- 
pressed in terms of scalar flux functions 



SB ±. A ^ J T 

VA—fT~ = ZX V_L* 
Bq 

and 

Su± = z x Vj_$. 
The ERMHD equations are 



(A16) 



(A17) 



(A18) 



A = -2- 



(A10) 



Note, the A m for the fast and slow roots can be trivially 
derived from equation (|A5|) since 5B y = 0. 



< T. 

(|Kingsep et al.1 Il990t 



2. Electron MHD 

Electron MHD (EMHD) is v alid for T, 
(ft ^> 1) and scales Aft > 1 



ISchekochihin et al"1l2009l ) and is most useful for describ- 
ing the whistler mode. At these scales, the ions decouple 
from the magnetic field and are treated as being station- 
ary within the framework of EMHD. Therefore, Ohm's 
law reduces to E = —u e x B/c. Inserting this form of 
Ohm's law into Faraday's law together with J = —neu e , 
we obtain the EMHD equation 



dB 

~dt 



V x 



(V x B) x 



cB 

Anne 



0. 



(All) 



The linear dispersion relation for this equation is 

ui = ±v A k\\dik, (A12) 
and the eigenfunctions for this system are 



6B X = --±6B 
k± 



and 



8B„ 



±l^5Bn. 

k_i 



Thus, the A m for the whistler root is 
k 2 + k 2 



•Am — 



k\ 



k\ 



(A13) 



(A14) 



(A15) 



<9$ 
~dt 



VA 



2 + ft [1 + 4» 

SB 1= , 
B V 

Sn e 
n 



z ■ V (p?Vi*) 



TiJ PiVA 

2 $ 



Vft PiVA 



(A19) 

(A20) 
(A21) 
(A22) 



where £ ■ V* = d* /dz + (l/v A )i ■ (Vj_* x V_l*). Note, 
the anisotropy of the system together with V • B = 
implies SB X <C <5-B z ~ SB y , so we need only determine 
the eigenfunction for SB y . 

The linear dispersion relation for this system is 



uj = ±1 



1 + Te/Ti 



^ 2 + ft(l + ft/T 4 ) 
and the eigenfunction is 

k±pi^ = PiVA-fT- = 
Bq 



kj_Pik\\VA 



(A23) 




1 + P 



(A24) 



<F 



Combining equations (|A20[) and (|A24[) yields the A r , 
for the kinetic Alfven wave. 



Am — 



2 + ft (1 + T e /Ti) 
ft(l + T e /Ti) ' 



(A25) 
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